Solving Block Low-Rank Matrix Eigenvalue Problems
نویسندگان
چکیده
To solve large-scale matrix eigenvalue problems (EVPs), a two-step tridiagonalization method using the block Householder transformation (HT) is often employed. Although based on dense arithmetic requires memory storage of O(N2) and an operations O(N3), in this study, these were reduced by approximating low-rank (BLR-) matrices. A special HT for BLR-matrices it are proposed to EVP with real symmetric BLR-matrix. In HT, vectors also formed BLR-matrices. It demonstrated how size m BLR-matrix should be determined confirmed that complexities O(N5/3) O(N7/3), respectively, typical cases when appropriate ∝ N1/3. numerical experiments string free vibration problem known analytical solutions, large eigenvalues, calculated eigenvalues converge toward ones accordance theoretical convergence curves. Owing complexity, was solved about N =300,000, which significantly larger than limit conventional matrices, within reasonable amount time CPU core. For calculation time, faster few tens thousands.
منابع مشابه
Stochastic algorithms for solving structured low-rank matrix approximation problems
In this paper, we investigate the complexity of the numerical construction of the Hankel structured low-rank approximation (HSLRA) problem, and develop a family of algorithms to solve this problem. Briefly, HSLRA is the problem of finding the closest (in some pre-defined norm) rank r approximation of a given Hankel matrix, which is also of Hankel structure. We demonstrate that finding optimal s...
متن کاملA Density Matrix-based Algorithm for Solving Eigenvalue Problems
A new numerical algorithm for solving the symmetric eigenvalue problem is presented. The technique deviates fundamentally from the traditional Krylov subspace iteration based techniques (Arnoldi and Lanczos algorithms) or other Davidson-Jacobi techniques, and takes its inspiration from the contour integration and density matrix representation in quantum mechanics. It will be shown that this new...
متن کاملSolving Overdetermined Eigenvalue Problems
We propose a new interpretation of the generalized overdetermined eigenvalue problem (A− λB)v ≈ 0 for two m × n (m > n) matrices A and B, its stability analysis, and an efficient algorithm for solving it. Usually, the matrix pencil {A− λB} does not have any rank deficient member. Therefore we aim to compute λ for which A − λB is as close as possible to rank deficient; i.e., we search for λ that...
متن کاملA geometric method for eigenvalue problems with low-rank perturbations
We consider the problem of finding the spectrum of an operator taking the form of a low-rank (rank one or two) non-normal perturbation of a well-understood operator, motivated by a number of problems of applied interest which take this form. We use the fact that the system is a low-rank perturbation of a solved problem, together with a simple idea of classical differential geometry (the envelop...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
ژورنال
عنوان ژورنال: Journal of information processing
سال: 2022
ISSN: ['0387-6101']
DOI: https://doi.org/10.2197/ipsjjip.30.538